Quantum interference of electromagnetic fields from remote quantum memories 
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We observe quantum, Hong-Ou-Mandel, interference of fields produced by two remote atomic 
memories. High-visibility interference is obtained by utilizing the finite atomic memory time in 
four-photon delayed coincidence measurements. Interference of fields from remote atomic memories 
is a crucial element in protocols for scalable generation of multi-node remote qubit entanglement. 
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Proposed approaches to scalable quantum information 
networks and distributed quantum computing involve lin- 
ear optical elements and single-photon detectors 0, 0, • 
Photoelectric detection events signal entanglement cre- 
ation and, by postselection, eliminate undesirable com- 
ponents of the electromagnetic field. While postselection 
has a residual negative effect on the scaling of the overall 
efficiency of quantum information protocols, this can be 
offset by quantum memory, a resource which provides the 
capability to perform quantum state transfer from mat- 
ter to light and vice versa, as demonstrated with cold 
atomic ensembles 0, 0, @ . These also act as sources of 
entangled photon pairs, with quantum memory enabling 
user-controlled delays between the photons. 

In order to distribute entanglement over a network con- 
figuration we must connect entangled elements at remote 
sites. This may be achieved by interfering photons, pro- 
duced at these sites, on a beamsplitter followed by co- 
incident photoelectric detection. The anticorrelation of 
coincidence counts is the signature of Hong-Ou-Mandel 
interference (HOM), whereby single photons simultane- 
ously incident at two input ports of a beamsplitter both 
exit in one or other of the output ports 0, 0- For 
distinct, remote quantum memory elements, HOM is a 
possible method for entanglement connection operations 
that scale efficiently with the number of elements. Sev- 
eral remarkable demonstrations of HOM using paramet- 
ric down- conversion (PDC) have been reported (see Refs. 
0, EJ and references therein). It has also been ob- 
served using photon pairs generated locally by a single 
source - a quantum dot an atom 0, and an atomic 
ensemble 13]. Recently HOM has been demonstrated 
with two (a) neutral atoms Q and (b) ions , in each 
case separated by a few microns. 

In this Letter we report HOM from remote sources 
(cold atomic ensembles) with quantum memory, located 
in adjacent laboratories and separated by 5.5 m (Fig. 1). 
In a particular ensemble signal photons are generated 
by Raman scattering of a write laser pulse with temporal 
profile if (t) (normalized to unity J dt \cp (t)\ 2 = 1), whose 
length is much greater than the ensemble dimensions. For 
an unpolarized ensemble of N atoms interacting with an 
off-resonant vertically (ey = — z) polarized write field 
propagating in the ^/-direction, the interaction picture 
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FIG. 1: Schematic showing Raman scattering of write pulses 
(W) at sites A and B with signal fields collected by polarizers 
Pi and P2 and optical fiber beamsplitter and directed towards 
detectors Dl, D2. A half- wave plate (A/2) may be inserted 
at site B to rotate light polarization. Raman scattering of 
delayed read pulses produces idler fields detected at D3, D4. 
The inset shows the atomic level structure and the write- and 
read- induced Raman processes. 



Hamiltonian for the system includes a term representing 
Raman scattering into the detected s igna l mode defined 
by an optical fiber, and described by 

H (t) = ih X ^ (t) (cos rj ft H (t) 4 + sin rj ft v (t) sty + h.c.(l) 

The annihilation operator for the Raman scattered 
field emitted from the ensemble is given by (t), 
A = H,V. These field operators obey the usual free- 
field, narrow bandwidth bosonic commutation relations 

?p\ (t) , ip{, (£') = 5 X ,\>S (t - t). The emission of H- or 
V-polarized signal photons creates correlated atomic spin 
wave excitations with annihilation operators sh,v given 
by, sy = —so, and 

sh = cos0s_i — sin0s + i, (2) 
where cos 2 = E m ^ X m-i/E 
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the spherical vector components of the spin wave are 
given by 
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An atom is prepared in the state |a,m) with probability 
Pm = 1/ (2F a + 1) for an unpolarized ensemble. Here 
| a) and \b) are the ground hyperfine levels and \c) is 
the excited level associated with the Di-line involved in 
the Raman process, with total angular momenta F a ^Fb 
and F c , respectively, and X m , a = C Fa 1Fc C^ b _ 1 a F ^ m is 
a product of Clebsch-Gordan coefficients. The spin wave 
is defined in terms of the /i-th atom transition operators 
a a m . b m , (t) and the write (j) w (r) and signal (j) s (r) mode 
spatial profiles (normalized to unity in their respective 
transverse planes) 



multi-mode, multi-particle system reduces to a discrete 
mode parametric interaction. 

Since d\ commutes with the Rayleigh scattering and 
undetected Raman scattering Hamiltonians, the Heisen- 
berg picture solutions are given by 
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where k s and k w are the detected signal and write beam 
wavevectors, respectively, and A is the effective overlap 
area of the write beam and the detected signal mode 

mm. 

The dimensionless parametric coupling constant x 18 
given by 



2d c bd ca \/k s k w n w N 7T^> ' 
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where A = ck w — (u c — uo a ) is the write laser detuning 
from the c <-> a transition, and d ca are reduced dipole 
matrix elements for the c <-» b and c <-> a transitions, n w 
is the average number of photons in the write pulse, and 
the parametric mixing angle rj is given by LLij 
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The interaction picture Hamiltonian also includes terms 
representing Rayleigh scattering, as well as the Raman 
scattering into undetected modes. One can show however 
that these commute with the signal Hamiltonian to order 
O and that they commute with ?p\ and s a . As 

a result, the interaction picture density operator for the 
signal-spin wave system (tracing over uncollected modes) 
is given by p = UpoW, where po is the initial state of the 
unpolarized ensemble and vacuum electromagnetic field 
and the unitary operator U is given by 

U = exp ^x cos? 7^L^L + X^ n V^v^v ~ h.c^j (6) 

with the discrete signal mode bosonic operator 
d\ =J dt if* (t) 7p\ (t), where A = H or V. In effect the 



= cosh (x cos rj) + sinh (x cos v) $h 
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These solutions allow calculation of the photoelectric 
detection signal for an atomic quantum memory ele- 
ment. In this work, two magneto-optical traps (MOTs) 
of 85 Rb, A and B, located in adjacent laboratories, 
serve as the basis for distant quantum memories (Fig.l). 
The corresponding Raman scattered signal field oper- 
ator for the detected mode with spatiotemporal mode 
<PA[B](t ~ z A[B]/c) is given by 
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with a similar expression for ensemble B. 

The fields from A and B are combined on a beam- 
splitter R + T = 1, where R and T are its reflectance 
and transmittance, and the fields J^a,i, in the out- 
put ports 1 and 2 are incident on detectors Dl and D2, 
respectively. We employ vertically (V) polarized write 
beams, derived from the same laser, and detect the hor- 
izontally (H) polarized signal fields, which are passed 
through polarizing cubes prior to the beamsplitter. 

The corresponding cross-correlation function 
G[, 12) (M + r) = (^iW^Hrj^lHr)^,^)) 
exhibits the HOM effect: 

Gf, 12) (M + r) = 



(7) 
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where s A[B] 



sinh(xA[£]COS77), and £ A[B] = 



y/hck a /(2€o)\^ 8l A[B] ( r ) I- Tne nrst , HOM, term on the 
right-hand side of Eq.Q exhibits two photon interfer- 
ence and can be understood in terms of conventional 
single-photon interference conditioned on the first pho- 
toelectric detection at time t 12]. For zero delay r = 
and a symmetric beamsplitter R = T = 1/2, this term 
gives zero contribution even for cp A ^ ip B . Alterna- 
tively, for if A = ip B it vanishes for arbitrary r. However, 

( 12) 

Gy (t,t + r) does not vanish completely due to contri- 
butions from multiphoton signal excitations (second term 
in Eq. Q)). To quantify the degree of the HOM effect, 
the following benchmark measurement is performed. We 
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insert a half-wave plate into the path of the signal field 
from ensemble B, rotating its polarization from H to V, 
thus nullifying the HOM effect. Quantitatively, in this 

( 1 2 ) 

case the corresponding correlation function G\_ ' (£, t-\-r) 
is given by Eq. J7J) , but now without the interference con- 
tributions (proportional to the product RT) in the HOM 
term. 

In our experiment magneto-optical traps (MOTs) of 
85 Rb are used to provide optically thick atomic ensem- 
bles at sites A and B (Fig.l). The ground levels {|a); \b)} 
correspond to the 5Si/2,-F = {3,2} levels of 85 Rb, and 
the excited level \c) represents the {bPi^^F = 3} level 
of the D\ line at 795 nm. For a linearly polarized write 
field we observe that the signal field is nearly orthogo- 
nally polarized, consistent with the theoretical value of 
cos 2 // = 91/122, Eq.©. 

In order to generate indistinguishable signal wavepack- 
ets from the two atomic memories, we produce their re- 
spective write fields by splitting a single pulse and direct- 
ing the outputs into identical 100 m long optical fibers. 
The two Raman-scattered signal fields produced at A and 
B are passed through polarizing cubes to select the H- 
components and coupled into the ends of a fiber-based 
beam splitter. The outputs of the latter are connected 
to single-photon counting modules Dl and D2. A half- 
wave plate is inserted into the path of signal field B (Fig. 
1) which allows us to vary the relative (linear) polar- 
ization of the detected fields. This allows us to detect 
parallel polarizations (||), which exhibit the HOM effect, 
and orthogonal polarizations (_L), which do not. 

Particular care is taken to eliminate possible sources 
of spectral broadening. Magnetic trapping fields are 
switched off after atomic collection and cooling, and 
the residual ambient field is compensated by sets of 
Helmholtz coils. All trapping and cooling light fields 
are switched off during data acquisition. The trapping 
light is shut off about 10 [is before the repumping light, 
preparing unpolarized atoms in level \a). 

In Fig. 2 we show the measured ratio of the photo- 
electric coincidence rates 7Z\\/7Z±, which are integrated 
over the duration of the write pulses. Our measure- 
ments exclude Rayleigh scattering on the write tran- 
sition by means of frequency filtering. The experi- 
mental ratio 1Z\\/1Z±_ is compared to the ratio of inte- 
grated correlation functions J J dtdrG^ (£, t + r) and 

J fdtdrG { l 2 \t,t + r), assuming identical wavepackets 
¥a = We observe scatter in the data beyond the de- 
viations due to photoelectron counting statistics. These 
indicate the level of systematic drifts encountered over 
several hours of data acquisition. 

The photoelectric coincidences arise from the signal 
field excitation pairs produced (I) one excitation from 
each ensemble; (II) both excitations from ensemble A; 
(III) both excitations from ensemble B. The HOM visi- 
bility of V = 1 — H\\/H± = 1/3 reflects the deleterious 
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FIG. 2: Ratio of measured two- fold coincidence rates for the 
J_ and || configurations). The parameter pi = (Ni + A^/iVr 
(averaged over the J_ and || cases). Here N±(N2) is the num- 
ber of photoelectric detections in detector D1(D2), Nt is the 
number of experimental trials. Theoretically it can be ex- 
pressed aspi = 6^4+^4) where ea(£b) ~ 0.05 — 0.07 is the 
overall probability to detect a signal photon from site A (site 
B) by either Dl or D2. Scatter beyond the estimated Poisso- 
nian level of uncertainty is consistent with systematic drifts in 
experimental conditions, in particular the single count rates 
from each ensemble. The solid line is our theoretical predic- 
tion based on Eq.Q, for R — T — 1/2 and cas\ = cbs 2 b . 

effects of contributions (II) and (III). These are relatively 
large because in the limit of weak excitation the spin 
wave-signal state is dominated by the vacuum contribu- 
tion. By detecting the presence of a spin wave atomic 
excitation in each ensemble, these contributions could be 
substantially suppressed, and the HOM visibility V — > 1 
in the limit that the excitation probability pi —> 0. 

We obtain high-visibility HOM fringes by means of 
a four-photon delayed coincidence detection procedure. 
This involves conversion of the spin wave excitation to 
an idler field by means of an incident read laser pulse 
which follows the write pulse by a programmable time 
delay St in the off-axis geometry |l8j : 6t is limited by 
the quantum memory coherence time r c [l^. By careful 
minimization of ambient magnetic fields, r c > 30 fis have 
been reported |20(. In this work we choose St — 100 ns 
in order to maximize the repetition rate of the protocol. 
The four-fold detection of the two idler and two signal 
fields involves HOM of the two signal fields and delayed 
coincidence detection of the idler fields at detectors D3 
and D4, as shown in Fig. 1. 

The four-fold coincidence rate is thus given by 

Tlf ~ s\sl{{R-Tf{l + 2s\){l + 2sl) 

+ 2RT (3s\ + 34 + 24 + 24)}- (8) 

We have again assumed identical wavepacket modes for 
both ensembles. 
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FIG. 3: Integrated four- fold coincidence rates /Wj 4 ^ and 

/W^ as a function of pi. Experiment, dots, theory, solid 
line, assuming identical signal mode wavepackets from each 
ensemble. Uncertainties are based on the statistics of the 
photon counting events. 

By inserting a half- wave plate into the path of the sig- 
nal field from ensemble B as before (rotating polarization 
from H to V), we suppress the HOM interference contri- 
butions, such that the four-fold coincidence rate becomes 

~ 44{(^ 2 + t 2 )(i + 24)(i + 24) 

+ 2KT(3s\ + 3s% + 2s 2 A + 2s%)}. (9) 

In separate sets of measurements we recorded pho- 
toelectric events with one, or other, of the two MOTs 
blocked, which allow us to determine the expected level 
of four-fold coincidences for orthogonal polarizations of 
the two signal fields (i.e., in the absence of HOM). 

In Fig. 3 we plot ft[ 4) /W^ and 71^ /w[ 4) along with 
the corresponding theoretical predictions. HOM inter- 
ference is manifested in that 7Z^ /W^ — > as pi — > 0. 

The highest observed visibility V = 1 - /w| 4) « 
0.86 ± 0.03. As the theory and the experimental data 
agree within the statistical uncertainties, this indicates 
very good wave-packet overlap of the signals produced 
by the remote ensembles. 

In conclusion, we have demonstrated quantum inter- 
ference of electromagnetic fields emitted by remote quan- 



tum memory elements separated by 5.5 m. Such high- 
visibility interference is an important element towards 
scalable distributed entanglement, the central resource 
in proposed quantum network and distributed quantum 
computing systems 
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